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Abstract 

The semi-classical collapse, including lowest order back-reaction, of a thin shell 
of self-gravitating quantized matter is illustrated. The conditions for which self- 
gravitating matter forms a thin shell are first discussed and an effective Lagrangian 
for such matter is obtained. The matter-gravity system is then quantized, the semi- 
classical limit for gravitation is taken and the method of adiabatic invariants is applied 
to the resulting time dependent matter Hamiltonian. The governing equations are in¬ 
tegrated numerically, for suitable initial conditions, in order to illustrate the effect of 
back-reaction, due to the creation of matter, in slowing down the collapse near the 
horizon. 


1 Introduction 


Much effort has been dedicated to the classical dynamics of collapsing gravitational bodies 
(see e.g. Refs. iii for early works and |Q] for a summary of more recent developments). 
In particular, the collapse of an isotropic homogeneous sphere of dust has been studied 
using the Arnowitt, Deser and Misner (ADM) construction Q both on relating the dust 
associated scalar field to time |P and including the quantum nature of the collapsing 
matter through a Born-Oppenheimer approach to the Wheeler-DeWitt (WDW) equation 
describing the coupled matter-gravity system (see i I, I and Refs, therein). Within the 
latter approach the evolution of quantum matter (dust) was studied both in the adiabatic 
approximation []R| and on relaxing it |]T]. 

Homogeneous dust in flat space-time without boundaries corresponds to a fluid with 
(possibly time-dependent) density and zero pressure and can be related to the homoge¬ 
neous mode of a free massive real scalar field. When dust self-gravitates and is confined 
into a sphere of finite volume, the metric inside the sphere is a section of a Friedmann 
universe (in the language of cosmological models) with constant spatial curvature and the 
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radius of the sphere depends on time Q. As we have previously indicated ©0 , us¬ 
ing the homogeneous mode of a scalar field to describe a sphere of dust is still plausible as 
long as one considers spheres whose radius is much greater than the scalar field Compton 
wave-length. In fact for such a case edge effects are negligible [^] . 

We have previously, through the use of adiabatic invariants [^] for time dependent 
Hamiltonians, relaxed the adiabatic approximation and thus allowed for the creation of 
matter (dust particles). However one will then eventually obtain a non-zero pressure and 
have deviations from the dust approximation. 

Rather than study an isotropic, but inhomogeneous, system it is easier to address 
the collapse of a gravitating matter shell which we now briefly review, also in order to 
introduce our notation. 

In general the dynamics of a system consisting of self-gravitating matter in a space- 
time manifold of 4-volume V with boundary B can be derived from the Einstein-Hilbert 
action with a matter term 


S= d'^x \/—g 

JV 


7^(4) 
16 vr G 


+ C 


M 


+ 


SttG 


d^xy/^lC= Ldr 


( 1 . 1 ) 


where G is Newton’s constant, g the determinant of the 4-dimensional metric, the 
Ricci scalar. Cm the Lagrangian density for matter and /C the trace of the extrinsic cur¬ 
vature of the boundary B on which a 3-dimensional metric with determinant g is induced. 

When matter is localized on an isotropic shell, we can introduce an adapted radial 
coordinate r such that the shell is parametrized by Xin < r < rout- Then the external 
geometry is uniquely determined by Birkhoff’s theorem to be given by ||l^ 


ds 


2 

in [out) 


-Ain{out) dt^ + dr"^ + dQ^ 


( 1 . 2 ) 


where = 1 for r < Xin, corresponding to a Minkowski metric inside the shell, and 
Aout = I — 2GM/x for r > Xoutj corresponding to a Schwarzschild metric of ADM mass 
M outside the shell. 

If the Schwarzschild radius of the shell, Rh = 2G M, is much greater than the thickness 
d ~ Xout — fin of the shell (as we shall show in the next section is true in our case) 
one can consider the “thin limit” d/xin —> 0 (with fixed matter proper energy), define 
R = {xin -|- Xout)(2- and the Lagrangian assumes the form |^, 


L = - 


RR , 1 /R 
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(1.3) 


where Pin (out) = \JAin (out) + R^, [a]lut = o(Dn) - a(w) for any function a(r) and a dot 
denotes the derivative with respect to the proper time r (dr = dt). The quantity m 
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is the volume integral of the matter Hamiltonian density and is equal to the proper energy 
of matter. 

The next step is the construction of the Hamiltonian, 


R 

H = P^R - L = -— {f3in - Pout) + m 


= m — 


V2R 




1 1/2 


(1.4) 


with 


Ph=^ = -- 

dR a 


ta.h- I 


(1.5) 


out 


It is a general feature of gravitational minisuperspace systems that the equation of motion 
is obtained on imposing the constraint H = 0, which expresses the invariance of the action 
under a time reparametrization. In our case the equation of motion takes the form Q 




I 2GM 




Gm 


( 1 . 6 ) 


By squaring twice this expression becomes 

( 1 - 

In the next section we consider gravitation with a minimally coupled scalar held and 
study under what conditions one obtains a collapsing shell-like solution. Once this has 
been achieved, we shall employ the ADM formalism in the Born-Oppenheimer approach 
to the WDW equation as was done for the sphere 0,0 and examine the effect of back- 
reaction due to particle creation on the shell during the collapse (section |3). This of course 
is tantamount to relaxing the adiabatic approximation. Lastly in section ^ our results are 
summarized and discussed. 


■ 2 GM G'^rrP 

~ “R 4R2 “ ^ 


2 The shell structure 


The purpose of this section is to illustrate how one can describe a collapsing matter shell 
in terms of a scalar held. In particular we shall investigate the conditions under which 
it is possible to neglect the thickness of the shell and use the thin limit Hamiltonian in 


Eq. (M) . 


Let us begin by considering an exterior Schwarzschild metric as given in Eq. (1.2) with 
Aout- The action for a massive scalar held <I> on such a background is 


S = ^ J dtdrr"^ 


1 - 


2GM\ 
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where [i is the inverse Compton wave-length of the scalar particle described by the field. 
The associated Klein-Gordon equation satished by a spherically symmetric (s-wave) scalar 
held solution <h(r, t) = ip{r) corresponding to a collapse is 


_l(,2_2GAfr)A + ^2,.2_ 


dr 


dr 


r -2GM 


(/? = 0 


( 2 . 2 ) 


On taking r 3> Rh^ one obtains 


dr'^ 


2d , 2 2n 2GM . 2 2n 

+ 

r dr r 


ip = 0 , 


(2.3) 


which of course resembles the hydrogen atom case. 

The classical limit for the Coulomb case has been studied in detail |1^|. In particular 
for large principal quantum number n and angular momentum I = n — 1 one obtains 
a probability distribution strongly peaked on classical circular orbits. For the case of a 
collapse (n large and I = 0) the probability distribution resembles a spherically symmetric 
collection of degenerate ellipses. This is analogous to the harmonic oscillator where for 
large quantum numbers the position probability density similarly approaches that for a 
classical oscillator of the same energy jl^ . This naturally suggests, in analogy with the 
harmonic oscillator, that one may consider a coherent state, that is a wave packet that 
will describe the collapse without spreading, at least for a sufficiently large period of time. 
Unfortunately for the hydrogen atom such a construction has not yet been successfully 
done (see e.g. p^, [^). 

An alternative way is to construct a shell composed of a suitable number of (isotropic 
I = 0) particles, bound together by gravitational interaction, whose mean position will 
follow the desired classical collapse. In order to illustrate such an approach let us first 
consider the semi-classical (WKB) limit in Eq. ( p.3| ) by setting (p ~ _ Retaining terms 

to lowest order in h one has 


dS \ ^ fc2/,,2 , ,2\ I 2G'i\T^2/,,2 


dr J 


+ + 


(/r^ - 2uj^) = 0 . 


(2.4) 


On setting {dS/dr)‘^ = (<C one obtains 


_ GM {tip) ^ ^2 - M^) 

2 {hp) r 2 {hp) 


(2.5) 


which for M = {hp)l2 is just the Newtonian approximation to the equation of motion for 
the collapse of a self-gravitating shell made of one scalar particle in s-wave with mass hp 
and potential energy G {hp)‘^/2r (see Eq. O))- 

The r.h.s. of Eq. (|2.5| ) is the (non-relativistic) energy for the motion of a particle in a 
gravitational field. We can generalize this to the case of N microshells, each one made of 
one s-wave particle with the same rest mass hp, constituting a self-gravitating macroshell 
whose mean radius satisfies a collapse equation similar to that for a shell of total mass 
h pN. We shall see that the essential difference with respect to one shell is that the 
N microshells form a radially localized bound state corresponding to a self-gravitating 
macroshell of finite thickness. 
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Let us start by considering the two microshell case. Each microshell, with mass hij, 
and radial position has a potential energy which depends on the self-interaction plus 
the interaction with the inner microshell (if present). Hence the total potential energy can 
be written as 


V{ri,r2) = 


2ri 


2r2 


G{h^if 


' 9{ri - ra) 
ri 


^ 6>(r2 - ri) ' 

r2 


( 2 . 6 ) 


where 6 is the step function. Since we are interested in a thin macroshell we introduce 
the mean radius coordinate R = {ri + r2)/2 and the relative coordinate r = ri — r 2 with 
|r| R. In this approximation it is possible to derive for the potential energy the following 
expression 


V{R,r) = 


G{2h^)^ 

2R 


Gjhfif 


r 


(2.7) 


For realistic values of the total mass of the macroshell (which means that if we wish 
to consider physical quanta of mass we must go to a many particle description with 
more particles on each microshell) the bound state has a characteristic time of oscillation 
which is several orders of magnitude smaller than the time needed to have a significant 
variation of the centre of mass position. For this reason it is possible to factorize the 
mean radius motion and study the quantum bound state adiabatically with respect to the 
slowly varying R{t). An interesting fact worth noting is the (slow) increase of the coupling 
during the collapse (since R decreases) which implies that the average shell thickness will 
decrease. For R constant the quantum problem for the relative r degree of freedom is well 
known and can be solved, with the spectrum of the allowed energy excitations being given 
in terms of the zeros of the Airy function or of its derivative. It is also possible to check 
that a semi-classical (WKB) analysis gives a good estimate of the energy levels, even for 
low quantum numbers. 

A similar conclusion to the above can be obtained in a classical general relativistic 
framework by using the junction conditions Q and taking the Newtonian limit. For a 
system of two non-crossing shells one has 


/ dri 
V dt 
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dr2 

dt 
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2GMi 
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2GMi 

r2 


M? GMi G^mf 

l + ^ +-^ + 

ri 4r{ 

(M2-Mi)2 ^ G(Mi + M2) , G^^ml 

1 ~\~ 9 H” ~\~ 9 

r2 4 r2 


( 2 . 8 ) 


where mi, m 2 are the proper masses of the shells positioned at ri < r 2 . Mi and t are the 
Schwarzschild mass and time coordinate for the metric between the two shells while M 2 is 
the Schwarzschild mass for the metric outside the outer shell. Since we wish to analyze the 
interaction between the shells in the Newtonian limit we neglect the terms proportional to 
G^ in the above. On taking mi = m 2 = one may sum (to this order) the two junction 
equations, multiply the result by hfi/2 and include the two position independent terms in 
the total energy E. One then has, for ri < r 2 , 



G{Ti^Y 3G(7l|u)2 

2 ri 2 r2 


(2.9) 
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An analogous term is easily obtained for the case ri > r 2 and finally the potential in 
Eq. (|2.6|) is re-obtained. 

If one wishes to consider a collapse approaching the horizon it is convenient to use the 
proper time of one of the two shells in Eq. (|2.8| ). In fact it is easy to check that on the 
horizon the ratio of the proper times of the two shells is 1 up to corrections of order d/R 
which we require to be negligible (see Eq. ( |2.21 ) for an estimate of this quantity). This 
implies that one is also allowed to consider the proper time associated with the mean shell 
radius. 

Let us now examine the more interesting N microshell case, again in the Newtonian 
limit. One may consider each microshell at a radial position r* and define rij = Vi — rj. 
The full interaction potential can then be derived: 


N 


V{ri,...,rN) = -J2 


2=1 


2ri 


-E 




2=1 


2 Vi 


E^(^b) 


( 2 . 10 ) 


where the first sum refers to the gravitational self-interaction and the second sum to the 
interparticle interaction. Introducing the mean shell radius R = one can write 

ri = R — {J2j^i^ji)/N and on using the thin shell approximation \rij \ <C R one has 


, G{Nhn)^ Gih^if^ 

E(ri,..., rjv) =-+ oEV E 


2R 


2R‘^ 

i<] 


( 2 . 11 ) 


In this expression for the potential energy the first term indicates that the mean shell 
radius motion follows approximately the Newtonian collapse of an object (the macroshell) 
of mass m = N hfi and the second term refers to the gravitational interaction of the internal 
degrees of freedom of the macroshell (relative positions of the microshells). One may 
study in the adiabatic approximation the quantum mechanical behaviour of the internal 
motion of the microshells but the problem is still complex and apparently unsolvable. It 
corresponds to the case of many particles with a confining potential linearly growing with 
the interparticle distances. This of course is analogous to the quark confining potential 
in hadronic physics, however for the latter case the total charge (colour) is zero. Further, 
due to the bosonic nature of the particles, in our case a condensate should be expected. 

At later stages of the collapse, as R{t) varies more rapidly, one expects non-adiabatic 
effects leading to transitions from the ground to higher (excited) states for some of the 
constituents of the particle condensate forming the macroshell. Clearly such excited states 
may decay emitting escaping radiation. However to describe this would mean introducing 
interaction terms and having an exterior Vaidya metric thus greatly complicating our 
simple model. The attitude we shall take and implement in the next sections is that 
once enough particles in the condensate are excited (thus widening the macroshell) they 


will collectively decay to the ground state by creating additional particles Thus the 
macroshell during its collapse will be globally described as a condensate of particles (all 
in the ground state) whose number varies (increases). 

In order to understand how the number N of microshells affects the global internal 
energy and the thickness of the macroshell we used a many body Hartree approach to eval¬ 
uate the energy for some trial parameter-dependent vacuum wave function. The average 
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Hamiltonian for a fully symmetric state is given by 


{H) = 




2 hfj, 


■ N dr 


d(p{r) 


+ 


dr 


N{N-l)G{nfif 

4^ 


GiNUiif 


J dr dr'\ip{r)\'^\r — r'\\(p{r')\'^ , (2-12) 


where we consider (p{r) = ce““ for the trial vacuum wave function with c 

a normalization factor and the phase factor describes the mean shell radius behaviour 
which is approximately free. Let us note that the asymptotic behaviour of the exact wave 
function for the N = 2 case is given by p = 3/2. One can verify that minimizing the 
average Hamiltonian for the relative motion with respect to a leads, independently of p, 
to the following qualitative behaviour 

{Hrel) -^npN{N , (2.13) 

and p only affects the magnitude of the proportionality constant. Further one also obtains 
an estimate of the width d of the packet 


1 fG{N-l) 
d ^ a ir r-^ — ' 






- 1/3 


(2.14) 


A similar behaviour can also be found if one calculates the energy levels semi-classically. 
Indeed for such a case the potential term may be evaluated by considering the microshells 
to be well ordered [23| so that the positions r^ = r^ — R (with respect to the mean shell 
radius) satisfy fi < f 2 < ... < tat and one will have > 0 with 


N N-l 


k=lj<k k=l 


N 

k{N -k) = 2 ^ ffc . 


k=l 


(2.15) 


The N microshell Hamiltonian then becomes 


N N 

« = E = E 


k=l 


k=l L 


\^^^rl + ^-^^^{2h)kk 


(2.16) 


We may now consider the fc-th microshell whose amplitude of oscillation will be the width 
d of the shell. Taking a closed cycle for the k-th. microshell one has a semiclassical quan¬ 
tization condition: 





- drk 



which leads to an energy 


Ek 



Gfip 
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and 
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e = J2e, 


k=l 


_/V5/3 

16 
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2/3 


(2.19) 


in agreement with Eq. ( p.l3 ). We note that the binding energy increases faster that the 
microshell number. In particular for a given n the smaller is R and/or the larger is the 
number N of microshells the smaller is the thickness d of the macroshell. 

Using natural units h = c = 1 one has G = ~ 10“^® GeV“^ and we can 

obtain some numerical estimates. If we consider matter quanta with a mass h/x = 1 
GeV and a number of quanta N = 10^^ (giving a total mass corresponding to 10“® solar 
masses), the Schwarzschild radius will be 2Gm = 10^®~ 0.2 cm and for a mean 
shell radius located at R = 10® Gm ~ 200 km, one has, from Eq. ( 2.14| ), an estimated 
thickness of the macroshell d lO^V^ ~ 10 ^ cm. With these numbers the potential 
energy of the mean radius of the macroshell is Gm?/{2R) ~ 10^® GeV, while the total 
internal energy for the other degrees of freedom, assumed to be in the ground state, 
is Eint ~ m{G N/ R^^/^ ~ 10®^ GeV (higher occupation numbers would lead to the 
previously shown n^^® behaviour). It is clear from these numbers that the condensate- 
maintaining mechanism we previously envisaged is plausible: a limited number of excited 
states can easily supply the energy for the creation of additional particles thus maintaining 
the condensate. 

From the above picture it is clear that, for a wide regime during the collapse, modelling 
the macroshell in terms of scalar fields confined around the mean shell radius is a good 
approximation. Thus we shall take for the scalar field describing one of the N particles a 
form <h(r,t) = (/>(t) (/?(r)/y7^ with ip{r) describing a packet with a width of order d 

around R. This time (p{r) will be normalized with the measure dr and not dr since we 
shall consider the dynamics of the field 0(t) with R{t) a time varying parameter. Then, if 
one considers a radial wave packet similar to the one used in the Hartree approximation 
and performs the radial integration in the action (O) for the scalar field, one obtains an 
effective action for i?i(t). In terms of the proper time r such effective action is given by 






fGNy/^ 


( 2 . 20 ) 


The last term in the integrand, which is related to the shell thickness, on comparing with 
the second term, can be safely neglected. Indeed, with the numbers previously given it is 
suppressed by a factor of 10^®. In the next section we shall use the formalism of second 
quantization in order to describe the N particle system discussed above. 

Lastly we note that, since we shall be interested in the evolution of the shell until it 
reaches the Schwarzschild radius Rh, which can be approximated extremely well by 2 Gm, 
we must check that the ratio 


d Eint (GN\^/^ 

R m \ R^ ) 

is still negligible with respect to 1 near Rh, that is 


( 2 . 21 ) 


G N ^ ^ Rh ^ 1 //i . 


For the values given above (/i = 1 GeV), we need iV 1/G/x^ ~ 10®®. 
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3 The quantized shell 


We may now write the WDW equation for the collapsing shell and the scalar field. The 
corresponding Hamiltonian will be given by Eq. (1.4) with m the average value for N 


particles of the Hamiltonian for the quantum scalar field obtained from Eq. ( 2.20 ). Let us 
now illustrate this. On expanding the Hamiltonian in Eq. (1.4) in G and retaining terms 
up to order G one has 


H = 




i?2 


+ R 


2 ±2 


p2 

- M H-^ 


2M 


GM^ 

2R 


= Hm + Hg , 


(3.1) 


where = R'^^/h and Pr = M R/vl + lP. Canonical quantization then leads to the 
WDW equation 


^2 . ,2^ ^2 52 GM ^ 

~ Jm ^ “ 2R 


and there is no operator ordering ambiguity. 


T(i?» = 0 , 


(3.2) 


3.1 Born-Oppenheimer approach 

We express the function 'I' in the form ^(R,4>) = ^{R)x{(l>,R) which allows us to obtain 
the WDW equation for the gravitational part 0] 


/ ^ GM2 92 

- M + 77TT + 


2R 




X, 


-0 


Hg + {Hm) ^ = 
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-ix 


dR‘^' 




92 


X, 


2M (x I x) 

where we have defined a scalar product 

(x I x) = Jd(f>x*i(p,R)xi(l>,R) 


2M\ dR^ 




and set 


with 






Further, we also obtain the equation for the matter field 

7 IH / £7 \1 ' , f ^X^ 7 
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5i?2 


92 

9i?2 


X ■ 


(3.3) 


(3.4) 


(3.5) 


(3.6) 


(3.7) 
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As yet there is no time parameter, however, on neglecting the r.h.s. of Eq. (3.3), one 
may introduce a time variable by taking the semi-classical (WKB) approximation for the 
gravitational wave function -0, 


'll) ~ 


dS, 


eff 


dR 


- 1/2 


gfi ‘^e// 


where the effective action Sgff satisfies the Hamilton-Jacobi equation 

( Hm ) + Hq = 0 


(3.8) 


(3.9) 


associated with the l.h.s. of Eq. (^.3|), and in the non-relativistic limit, .R <C 1, it is given 
by 


Seff = l f dr 


2M + M + 


GM^ 


R 


-2{H 


M . 


(3.10) 


where (Hm) is evaluated on the semi-classical trajectory where -0 has support. 
The explicit expression of Eq. 


R2 = 


M 




+ 


GM 


(3.11) 


may be compared with the classical Eq. (O- The two equations lead to the same evolution 
if 


( Hm)= m 


(3.12) 


and 


G^' 


<C 


4R2 " R 

2 

—^ — 1 ~ — (M — m) 


Gm? 

RM 


< 1 


m/M ~ 1 . 


M 

With the above conditions we can obtain as our equation of motion for R 


R^ - 


GM m 2 


- 1 


(3.13) 


(3.14) 


which is just the equation for radial geodesics in Schwarzschild space-time. In particular, 
on comparing Eq. ( |3.14 ) with the standard notation for a radial trajectory with geodesic 
energy E in a Schwarzschild space-time with mass parameter M (see e.g. |12|), 




(3.15) 


we find M = M/2 and E = M/m. Hence the shell evolves as a particle in a Schwarzschild 
manifold of total energy M/2, that is half of the total energy of the shell. The solu¬ 
tions to Eq. dUl) are known. Here we only recall for later reference the case m = M, 
corresponding to the separatrix between bound orbits and scattering trajectories: 


Rc{t) = Ro 


1 - 


IGM 

14 


2/3 


(3.16) 
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where Rq = R{0). 

The proper time variable associated to the WKB gravitational wave function is given 
by the relation 


d'ljj d 

i dSeff 1 

( dSeff \ 

1 

to 

dR dR 

n dR 2 

\ dR J 

dR^ 


i.± 

^ dR 


Z - . 7 9 1 

=- M lb — - 

h ^ dr 2 


dS. 


eff 


dR 


-1 


d‘^s, 


A 


(3.17) 


If the r.h.s. of Eq. (3.7) and the second term in Eq. ( 3.17| ) are negligible, one obtains a 
Schrodinger equation for the matter wave-function Xs, 


i h 


dXs 

dr 


2 






+ R 


■2 i ,2 


= Hm Xs 


(3.18) 


where we have scaled the dynamical phase 


Xs = X exp 


i 

h 


{Hm{t') ) dr' 


(3.19) 


3.2 Adiabatic invariants 

In order to study quantum systems with explicitly time-dependent Hamiltonians we can 
use the method of invariants [14|. Given a Hamiltonian a Hermitian operator /(r) 

is an invariant if it satisfies 


• ^ dI{T) 

OT 




= 0 . 


The general solution Xs('^) to the Schrodinger equation 

■ ^dxs{T) fr r \ r \ 
i Tb = Hm[t) Xs(t) 

can then be written in the form 

I X, r)/s = I n,T)j , 


(3.20) 


(3.21) 


(3.22) 


where | n,T)j is an eigenvector of /(r) with time-independent eigenvalue An and Cn are 
complex coefficients. The phase ^niT) is given by the sum of the geometrical and dynam¬ 


ical phases associated respectively to Eqs. (|3.6|) and (|3.19|) , 

7 /-T ^ ^ 

Vn{r) = - J I {n,T' \ hdr' + i Hm{t') \ n, t' )j dr' . 


(3.23) 


The Hamiltonian in Eq. ( p.l8 ) describes an harmonic oscillator of hxed frequency ^ 
and variable mass R? /ji. Thus we can introduce the linear (non-Hermitian) invariant 


4(r) = S(r) 


(3.24) 
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with the phase 0(r) given by 


0(r) = f dr' 
J rn 


R? (r') 


(3.25) 


^ I 

b[T) = .— —h ^ \ xns- X — 
X \ ^ // 


The function x(r) is a solution of 

■■ n ■ R 2 

^. + 2 :.- + ,. x = 
with suitable initial conditions. 

The system admits an invariant ground state (vacuum) | 0, defined by 

4(t) I 0,r)f, = 0 , 

from which one can define a basis of eigenstates .6 = {| n, r )^}, 


(3.26) 


(3.27) 


(3.28) 


\Jn\ 

= ^ I 0,r), = I n,T), , 

Vra! 

where (pQ is arbitrary and can be replaced by 0/2, and 

h\ n,T)f, = y/n\ n- l,r)f, 
b^\n,T)^ = Vn + 1 1 n + 1,r)^ . 


(3.29) 


(3.30) 


Since b,b^ =1 we will refer to b^ and b as the invariant creation and annihilation 
operators. We can then introduce the Hermitian quadratic invariant operator 


R^filR + ^\=n(pb + ^ 


(3.31) 


in which b^ b is the invariant number operator. 

We also define the particle creation and annihilation operators and d (with d 
1) as 


a(r) = 
a'''(T) = 


I • 


Ry ’ 


(3.32) 



in terms of which the matter Hamiltonian operator can be expressed as 


Hm — h/i ( + — 


(3.33) 


where iV = a is the particle number operator which counts the number of quanta of the 
scalar field (j). The corresponding vacuum state | 0, r)^ is defined by 


a| 0,r)„ = 0 

and a complete set of eigenstates can be given as ^ = {| n, r)^}, 

(at)' 




0>^)a t 


with 


«l n,T)^ = ^1 n - l,r)^ 


a' n 


,T)a = VnTT\ n + l,r)„ . 


The two Fock basis A and B are related. On using 


= vt® (^+ 


"^0 


we obtain 


where 


f - * 

a = B*b + A*b^ 

-V /s or 

= Bb^ + Ab 


Mr) = -R [x- 
B{r) = ^R (x + 


-b]+xf 


(b^ + b 


b = B a — A* 
b^ = B*a^ -Aa 


1 . X 

— i 


R^ X /i 

1 X 

— i 


(3.34) 


(3.35) 


(3.36) 


(3.37) 


(3.38) 


(3.39) 


R? X n, 

are the Bogoliubov coefficients. The two basis A and B will coincide for r = tq if 

6(ro) = a(ro) ^ /r/c(ro) = Hm(to) , (3.40) 

that is H(ro) = 0 and B{tq) = 1, which means that for the function x{t) we must require 


x(ro) = l/R 
x(ro) = 0 . 


(3.41) 


In the following we shall take tq = 0 so that the adiabatic approximation holds at the 
beginning of the collapse. 
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If we define a = Rx, Eq. 


takes the form 


— K + fl cr — —^ , 

which is the equation for a harmonic oscillator with time-dependent frequency fl, 

fl2 = i+ ^ 


;u2i?2 


On using the equation of motion ( |3.14| ) one has 

GM 


~ 1 


2/i2 R3 


= 1 + 


2 R^ ’ 


(3.42) 


(3.43) 


(3.44) 


where 6 (also suggested by Eq. (|3.16| )) plays the role of a parameter measuring non- 
adiabaticity since for <5 —> 0 the adiabatic limit is recovered. On expanding in 6, an 
approximate analytic solution to Eq. (3.42) is given by |24] 


8 


Then 


a 1 

= 7? - 7? 

It rC 


8 RS 


and 


0(r) ^ /i / dr' 

Jo 


1 + ^^ 
^ 4 i?3 


(3.45) 


(3.46) 


(3.47) 


The introduction of the eigenstates of adiabatic invariants allows us to examine the 
particle production induced by the evolution of the radius of the shell. As we have seen 
previously, in order to have a physical shell we must consider a large number of particles. 
Thus let us consider a linear combination of eigenstates of the invariant number operator 
defined in Eq. (|3.29|) , 




1 


N+X 

V2A._? 


(3.48) 


n=N-X 


with 1 <C A <C A^. Then we can evaluate the quantities of interest as expectation values 
for this state. On defining (O) = bs{J^x,T \0\ we obtain (up to corrections of 

0{\/N)) 


(<(') 

(^<A> 


V2hN X cos 0 
xR^ 


V2hN 


1 


cos 0-sin 0 

// X 


(3.49) 
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and 


2 


((Ac^) 

((Avr^)^ 


— — X 


2 I x2 ^ 


which lead to an uncertainty relation 


/ \2 1/2 

((a.^) ) {(Air*)") 


2 , 1/2 h L 2-2 

1 H- TT x^ . 


2 V 

On substituting the explicit expression for x and x from Eq. (3.46) we obtain 


((a<^)' 




h 1 


6^ R?' 


^ 2 R^ V A R^ 


h 

2 




i?3 


and 


/ ''\2 1/2 „ 
{(A0) ) ((Ari*)") 


2 ,1/2 h 

“ 2 


JL ^_^A2< 


1 + 


1/2 


We also find 


and the expectation value of the Hamiltonian is 


{ic)c^n {n + 


1 


{Hm ) — 


_h_ 

2i?2 




On approximating 


AA,r)(,^ ~ I /V,r) 


we have 


(ffM) =r ^(iV + i 


65 


1 o „o 


a;2 i?2 


+ i?2 + 




~ hfiN 


1 + 




i2 


= mo 


1 + 


2/r2 7^2 
i?2/23 


Rh i?2 
= mo + Am{T) , 


(3.50) 

(3.51) 

(3.52) 

(3.53) 

(3.54) 

(3.55) 

(3.56) 


(3.57) 


to leading order in 6. This value can be identified with a varying proper mass m^r) of the 
shell. 
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3.3 Consistency conditions and back-reaction 


The previous results were obtained on neglecting the r.h.s.s and the higher WKB order 
terms in the equations for the matter and gravitational wave-functions. Thus, in order to 
check the validity of our results, we must estimate the neglected terms and compare them 


to the corrections we obtain with respect to the classical trajectory in Eq. (3.It). In the 
following it will be sufficient to consider states | the approximation (3.56) and 

study their evolution for the mean shell radius R greater than the Schwarzschild radius 
Rh- 

The r.h.s. of the gravitational equation (|3.7|) can be estimated according to 


2M 


,{N\dl\N), = - 


4,,2 


TiV 


4Mi?2 


m 


E 

P^N 


{N |(at2 + a2)| P){P I (at2 +a2)| N) 
{En - EpY 




(3.58) 


where we used for the time-dependent number of scalar quanta the relation N{t) = m/hfi. 
Therefore, the above correction to the Hamilton-Jacobi equation for gravity is negligible 
with respect to the matter source m = ( Hm ) if 


R > l//r , 


(3.59) 


which means that the fluctuations associated with the quantum nature of matter be¬ 
come dominant for a radius of the shell of the order of the Compton wave-length of the 
constituent quanta (which appears as a general property of quantized models, see also 
Refs. [^, 1^). It is worth noting that the condition ( 2.22| ) which ensures that the width 
of the shell d <€. R coincides with (|3.59|) when R = Rh- 

For suitable initial conditions it is possible to have R > Rh 1//4, so that one may 
follow the evolution of the shell down to its Schwarzschild radius while using the semi- 
classical Hamilton-Jacobi equation. However, on using the second line in Eq. ( 3.571) to 
estimate Am one finds that the ratio p between Am and the quantity in Eq. (|3.58|) is 


p< 2M^ 


(3.60) 


where the maximum (of order 2) is obtained on the horizon. Hence, in order to compute 
consistently the back-reaction on the trajectory R{t) due to an increasing number of 
quanta, one must include the contribution in Eq. ( p.58 ) into the Hamilton-Jacobi equation. 
This leads to the equation of motion given below by the first expression in Eq. 


The r.h.s. of the matter equation ( |3.7| ) induces quantum transitions from the (otherwise 

L)^ with L ^ N. It is then useful to estimate the 


conserved) state 
projection 


N)^ to states 


2,,2 


b{L\dl\N), = 


R2 


E 

P^L,N 


{L\{a^^ + d^)\P){P\{d^^ + d^)\N) 
{Ep - El) {Ep - En) 


m 


Ap?R^ ’ 


(3.61) 
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for L = ± 4. Since the latter expression leads to a r.h.s. in Eq. (3.7) of the same 


order of magnitude as Eq. (|3.58|) , it also leads to the same condition, Eq. (|3.59D , for the 
validity of the Schrodinger equation ( p.l8 ). Further, in Eq. ( 3.18| ), if one wishes to include 
corrections related to Am, one must also include the contribution of the r.h.s. of Eq. (3.7) 
as was done above for the gravitational part. In the following we shall check that, for 
physically plausible data. Am <C mo and both the above corrections can be neglected. 

Finally we must check the consistency of the WKB approximation used in Eq. ( ^.17 ). 
For the matter equation (3.7) the correction coming from the (neglected) prefactor in the 
semi-classical wave function is approximatively given by 


Pr 


n^R 


< 


2M RP^ 2MiP ~ 


(3.62) 


where the last expression is the value on the horizon Rh- For the gravitational equation 
the corresponding correction is given by 


Jm 


3 P^R 


PrR 


+ 


Pr 


4i?2p2 2R^Pr 2R^Pr 




2M 


5R^ 1 

4 .R4 + 2 ^ 


< 




G2 M3 


. (3.63) 


Both corrections above, when compared with the matter source m are negligible provided 


Rr ^ h/M , 


(3.64) 


where h/M is the Compton wave-length of the shell as a whole and is much smaller than 
l//i (see Ref. [|T^]). Thus if condition (p.59D is satisfied, so is (|3.64|) . Further, the first 
correction is proportional to N~‘^ and the second one to N~^. From the expression for 


m in Eq. ( 3.57 ) one can infer that Am is proportional to the ratio between the Compton 
wave-length l//i and the radius at which Am is computed (in our case R > Rr)- Thus 
Am ~ 1/N (on the horizon where R is a constant) and both corrections are negligible 
with respect to Am for N ^ 1. 

In order to compute the back-reaction on the trajectory of the shell, one can inte¬ 
grate the system of coupled ordinary differential equations governing the semi-classical 
(Newtonian) trajectory, which we denote by Rg, and the corresponding proper mass m. 


d 2 _ GM I n f-I rn\ m? 

— Re. ^ ^ P M) TJPW^ 

(1+2)5^, 


Rs 
m = mo 


(3.65) 


The last term in the first equation above corresponds to the quantity in Eq. ( ^.58 ). As 
initial condition we shall set m(0) = mo = M. It is then natural to compare the result 
with Rc in Eq. ( 3.16| ) which corresponds to the Newtonian limit for the trajectory of a 
thin shell with m = mo = M constant. 

We have computed Rs{t) and m{T) numerically for a set of initial number of quanta 
N > 10^®, which is the limit required by our model of the shell described in section (see 
Table the third case is included only as a limiting example). The general picture is 
that Rs collapses slower than Rc, however the difference Rs — Rc computed at the proper 
time when R^ = Rr is a fraction of l//i and decreases for increasing N. Therefore, in our 
approach, the difference between the semi-classical trajectory Rg and Rc is negligible in 
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N 

Rh 

AN 

Rg - Rc 

Rs-Rc. 

Rc. 

rri^ 

8A?^ 

O 

r—1 

500-^ 

4 

• 10^4 

1 

O 

r—1 

CN 

4•10"^ 

0.13 

4 • lO^y 

50-^ 

A* 

4 

• 10^^ 

2 • 10”^ -f 

4 • 10-'^ 

0.13 

4 • 10^« 

5-1 

_ ^ _ 

4 

7io35“ 

2•IQ-i A 
_ 

4•10"^ 

0.13 


Table 1: Samples of numerical estimates of the effects associated with non-adiabaticity. 
N is the initial number of particles (with Compton wave-length Ijfj, = 2 - 10“^'^ cm) in the 
shell; Rh the corresponding Schwarzschild radius; AN the number of particles produced 
at the proper time at which the shell coordinate Rc reaches the horizon; for the same 
proper time the relative position of the shell when back-reaction is included is given by 
Rg — Rc] finally, in the last column we report the ratio between the relativistic correction 
and the Newtonian potential on the horizon, (see Eq. (3.13)). 


all physical situations. Further, relativistic corrections associated to the neglected term 

m?/A R1 in the equation of motion (O) amount to a deviation from Rc of about 13% 
on the horizon (see last column in Table |I]) which is much greater than the computed 
relative deviation of the trajectory due to a changing mass or to quantum fluctuations 
(fifth column in Table |I]) . Therefore one can safely approximate Rg with Rc down to the 
Schwarzschild radius. 

The production of matter instead is appreciable. In the second case displayed in 
Table ||, the number of produced particles AN = Am/hfj, is of order of 10“^ times the 
initial mass (see also Fig. |l]). In general the data shown in Table || are in agreement with 
the expected behaviour AN oc N~^. 

Lastly for Rg ^ 0 the number AN diverges and our approach breaks down completely 
since all the consistency conditions are violated. 

4 Conclusions 

Our purpose has been the description of the gravitational collapse of a self-gravitating 
shell of quantum matter. We have attempted to obtain a “realistic” description of such a 
shell by considering a bound state of a sufficiently large number of scalar particles. It was 
heuristically found that it was possible to obtain realistic and “thin” shells. 

Since our shell consists of scalar particles it is plausible to assume that the particles, at 
least initially, will all be in the ground state and form some sort of condensate. Subsequent 
evolution during the collapse will cause transitions to higher excited bound states which 
presumably, once a sufficient number of them is formed, will decay collectively through the 
creation of an additional particle in the ground state. This will lead to a condensate with 
an increasing number of particles. We have pointed out that this is somewhat analogous 
to the hadronization process in QCD. 

The above approach then led to an effective time dependent Lagrangian for matter. 
On including the shell Lagrangian and quantizing the matter-gravity system we obtained 
the corresponding Wheeler-DeWitt equation. On subsequently considering the semiclas- 
sical limit for gravity and neglecting fluctuations, we introduce a matter Fock space and 
reproduce the precedingly obtained Einstein equation for a many scalar particle shell to¬ 
gether with the Schrodinger equation for matter. The latter system was then solved using 
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Figure 1; Evolution of the radius Rs/Rh and relative mass 10^ (m — M)/M of the shell in 
proper time for the second case in Table 1 and Rs{0) = WRh- 


the method of adiabatic invariants. 

On analyzing our results, the main conclusion is that, for a suitable number of initial 
particles constituting the shell, the matter proper energy increases appreciably when the 
shell radius approaches the horizon. This corresponds to an increase of the number of 
particles, that is matter creation. Correspondingly, one obtains corrections to the equation 
of motion for the mean shell radius which lead to a slow down of the collapse. However, for 
the range of parameter considered, this change does not produce observable effects on the 
trajectory of the shell. It is clear that suitable conditions should lead to more significant 
effects, however our approach would then break down. 

It would be of interest, within our quantitative approach, to consider an extension to 
many shells leading to a more general matter distribution in order to study its collapse, 
again taking into account matter quantum effects. We hope to return to this point. 
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